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The Casimir effect results from alterations of the zero-point electromagnetic energy introduced by 
boundary-conditions. For ferromagnetic layers separated by vacuum (or a dielectric) such boundary- 
conditions are influenced by the magneto-optical Kerr effect. We will show that this gives rise to a 
long-range magnetic interaction and discuss the effect for two different configurations (magnetization 
parallel and perpendicular to the layers). Analytical expressions are derived for two models and 
compared to numerical calculations. Numerical calculations of the effect for Fe are also presented 
and the possibility of an experimental observation of the Casimir magnetic interaction is discussed. 



I. INTRODUCTION 

Since its discovery, the Casimir effect has gradually 
become a much-discussed subject in physics. Originally, 
one understood by the Casimir effect the attractive force 
between two metal plates in vacuum as a result of zero- 
point quantum fluctuations |Q. Nowadays the term is 
used for a much broader range of effects, all involving 
the influence of boundaries on fluctuations. As such, the 
Casimir effect plays a role in quantum fleld theory, atomic 
and molecular physics, condensed matter physics, gravi- 
tation, cosmology and so on. A thorough review of the 
Casimir effect in all these fields was published recently [|j . 

For two uniformly magnetized ferromagnetic plates 
held parallel to each other, it is shown in a previous pa- 
per 1^ that the interplay of the Casimir effect and the 
magneto-optical Kerr effect gives rise to a new long-range 
magnetic interaction. In this magnetic Casimir force 
was found to decay with interplate distance D as D^^ 
in the limit of long distances, and as for short dis- 
tances. In this case, the ferromagnetic plates were de- 
scribed with a Drude model and the magnetization was 
defined to be perpendicular to the plates. In view of 
future experimental investigations of this new magnetic 
Casimir force, it would be useful to study the case where 
the magnetization is parallel to the plates since this sit- 
uation is easier to obtain in an experimental setup. This 
subject will be st udied in the pres ent paper and a force 
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which decays as D^^ in the long-distance limit and as 
D^^ in the limit of short distances is found when the 
Drude model is used. This behavior is interesting since 
it means that the force is larger, and thus easier to mea- 
sure, for in-plane magnetization than for perpendicular 
magnetization at sufficiently small distances. Next to the 
Drude model, another more realistic model is also stud- 
ied. In this so called hybrid model, a plasma model is 
used for the diagonal element of the dielectric tensor of 
the magnetic plates, and a single absorption line model 
for the off-diagonal element. As for the long distance 
limit, the force in this model goes like for the case 
of perpendicular magnetization and as _D"^° when the 
magnetization is parallel to the plates while the behav- 
ior in the short-distance limit is unchanged. Finally we 
will present some numerical calculations of the interac- 
tion for Fe, in which experimental data for the elements 
of the dielectric tensor are used. An experimental setup 
to measure the magnetic Casimir force is also discussed. 



II. GENERAL THEORY 

Consider two uniformly magnetized ferromagnetic 
plates of infinite lateral extension held parallel to each 
other. The Casimir interaction energy per unit area at 
T = between the plates can be expressed as M : 



(2^)^ Jo 



+ 00 



dw / d^fcii Im Trln(l - R^Rse 



(1) 



where k±^ and fc|| are the components of the wavevector 
perpendicular and parallel to the mirrors. The 2x2 
matrices Ra and R^ contain the reflection coefficients of 



the two mirrors: 
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The index s (resp. p) corresponds to a polarization with 
the electric field perpendicular (resp. parallel) to the 
incidence plane. We will adopt here the usual conven- 
tion that the s axis remains unchanged upon reflection. 
Since the reflection coefficients are dependent on the di- 
rection of the magnetisation of the mirrors, it is clear 
from Eqs. (|l|) and (^) that the magnetic Casimir energy 
between the mirrors will differ for the situations in which 



the magnetizations of the two mirrors are parallel (FM) 
or antiparallel (AF). This will result in a net magnetic 
Casimir force per unit area AJ-' = Taf ~ m between 
the mirrors, different from the ordinary Casimir force dis- 
cussed in Ref. Q|. 

If a change of integration variables {lo, k\\) — > {uj, k±, Lp) 
is performed and complex integration methods are used 
as in [Q, Eq. ([|) can be written as: 
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In general, the reflection coefficients contain terms of dif- 
ferent orders of the magneto-optical constant Q. In our 
calculation, only terms up to first order in Q will be con- 
served. When the magnetisation direction is reversed, 
these terms will change sign. Since the first order terms 
are usually much smaller than 1 and than the terms that 
are independent of Q, it is possible to expand expres- 
sion (0) to lowest order in the linear terms. 



A. The Polar Configuration 



After some algebra, we find for the situation with mag- 
netisation directed perpendicular to the plates (we will 
call this the polar configuration from now on): 
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where the reflection coefficients have to be evaluated at imaginary perpendicular wavevector and frequency. In this 
equation, the reflection amplitudes are supposed to be identical for the two mirrors. Otherwise, the squared reflection 
coefficients have to be replaced by the product of the coefficients for the separate mirrors (e.g. r^^ ^i^£)- The 
integral over the angle ip is already performed. The reflection coefficients for a mirror in the polar configuration are 
given in as: 
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with: ^ = A/a;2(e^^(iu;) - 1) + (fc_Lc)2 
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B. The In-plane Configuration 

For the case where the magnetisation is parallel to the 
plates (we will refer to this situation as the in-plane con- 
figuration from now on), not only r^p, but also Vpp will 
contain a term that is linear in the magneto optical con- 
stant. As a consequence we find two contributions to the 



Casimir magnetic interaction energy. The first one (A£j ) 
results from the longitudinal Kerr effect, while the second 
term { A£\ ) is a consequence of the transversal Kerr ef- 
fect. Again the integral over ip can be performed directly. 
We obtain (for identical mirrors): 



J 



3 



Pi 


n + oc 




I 






47^2 


I 

Jo 


h 


r + oo 


~^ . 


la 


-h 


p + CO 


2^ 


Ja 



dk [ k \ 



dk±k± 



dfc_L kj 



k± c 



do; Re 



duj Re 



duj Re 



(1 - r2^e-2fc-Li3)(i _ J.2 e-2fc_LD) 







i'k±c 

t]_ / dtjRe 
Jo 



At-2 A~2fcj_-D 

(1 - r2pe-2'=^^)2 

(rjp)2[l-r2^r2^e"4fc^^]e-2fc^^ 
([1 - r2^e-2'=^^][l - r2pe-2fc^^])2 
Ar2p(l + r2pe^2fc,^)g-2fc,i3 



(1 - r2 e-2'=^^)3 



(6a) 
(6b) 
(6c) 
(6d) 



r 



Of course AfH = ASf + AfJ and AJFH = AJ^j' + AJ^|. integration over the angle </? is already performed. The 
For two different mirrors, the squares of the reflection reflection coefficients in Eq. (|^) are given in jsj: 
coefficients have to be replaced as mentioned above. The 
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with again: ^ = ^^2 (£3,3,(10;) - f ) + (fcj_c)2. rss{iijJ,ik±) 
and rpp{iuj, ik±) are still given by Eq. (|5a|). Note that the 
contributions arising from the longitudinal and transver- 
sal Kerr effect are of opposite sign and therefore tend to 
cancel each other. On the basis of Eqs. (^0), we will cal- 
culate the Casimir magnetic energies and forces for two 
simple models in the next two sections. In section 
we will use these equations to numerically calculate the 
interaction for iron plates. 



III. THE DRUDE MODEL 

Consider two identical magnetic mirrors with a dielec- 
tric tensor described by the Drude model: 
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In this equation Up is the plasma frequency defined by 
ujp = Anne'^/m*; ujc is the cyclotron frequency given by 
LOc = eBcs/rn*c where BcS is the effective magnetic field 
experienced by the conduction electrons as a result of 
exchange and spin-orbit interactions; r is the relaxation 
time. In the usual situation: lUcT ^ 1 <C ujpT. 



There are three important distance regimes to con- 
sider. In the long distance limit {D 3> cr) the dominant 
part in the integrals in Eqs. (||) and (||) comes from the 
region w < k±c « c/D <C l/r. In this range, one has: 

£xx{ii^) ~ £xx{^^) - f ~ > 1, (9a) 
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With these approximations, one finds for the reflection 
coefficients: 
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For the polar configuration, we arrive at: 
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FIG. 1: Absolute values of the magnetic Casimir force (per unit area) with the mirrors described by a Drude model. Numerical 
results (solid curves) are compared with the analytical expressions (dashed curves) for the polar configuration (a), the in-plane 
configuration (b), the term A.7^|' resulting from the longitudinal Kerr effect (c) and the term from the transversal Kerr effect 
A^jl (d). 



While for the in-plane configuration it is found that: 
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The second regime is that for intermediate distances 
(c/wp < £) < cr). Now the integrals in Eqs. (^) and (||) 
are dominated by the range 1/r <C a; < A;^c « c/D ^ 
ujp. For the elements of the dielectric tensor one then 
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In this case, the reflec tion coef ficients r, 
Arpp still satisfy Eqs. (10a- 10c| ), while 
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Since the reflection coefficients for the in-plane configura- 
tion in this regime are not different from the ones in the 
short distance limit, the expressions (^2|) are still valid for 
the Casimir magnetic energies and forces in the in-plane 
configuration. However, for the polar configuration, one 



5 



has: 



1 he LO^ 

24 



UJv 



1 he uj^ 

8 



(15a) 
(15b) 



The third regime to be considered is the hmit of short dis- 
tances (-D <C e/ujp). Here one has to distinguish between 
two regions: (i) uj < k±e ^ cOp and (ii) tOp 'ti uj < k±c. 
In region (i) the dielectric tensor elements are given in 

Eqs. ([l3|), so the reflection coefficients r^, rip and Arpp 
are the same as those in the intermediate distance regime, 
but we will now make an expansion for r^g and Vpp around 
— 1 and 1 respectively: 
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In region (ii) the dielectric tensor elements are given by: 
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With these approximations for the reflection coefficients 
in regions (i) and (ii), one finds the following expressions 
for the energies and forces (only the dominant term is 
given) : 



From these expressions, it is obvious that A£» « AfJ, 

and AJF" « ^-^J ^o^' distances small enough. Note that 
for the polar configuration, the exponent of the depen- 
dence with respect to D obtained here in the short dis- 
tance limit differs from the one obtained in Ref. |^ . This 
is due to the effect of multiple reflections, which were ne- 
glected in Ref. ||^ in this regime. It is interesting to note 
that although the reflection coefficients are much smaller 
than 1 in this high-frequency limit, the effect of multiple 
reflections is so important that the analytical dependence 
with D is modifled. This is a unique feature of the mag- 
netic Casimir effect. 

It is clear that in the polar configuration, the energies are 
always negative. For the in-plane configuration however, 
AE\ is positive, while AE\ is negative, so the sign of the 
resulting energy A^" = ASf + A£\ will depend on the 
magnitude of these two terms in the different regimes. As 
a result, a change of sign of the interaction is observed; 
in the long and intermediate distance regimes, the total 
energy Af" is negative, while for short distances it is 
positive. So whether the magnetic Casimir interaction is 
negative or positive depends on the distance between the 
mirrors. 

We numerically calculated Eqs. (||) and (||) and com- 
pared them to the analytical expressions derived above. 
Details of the numerical procedure will be given in sec- 
tion 0. The absolute values of the magnetic Casimir 
forces per unit area (both numerical and analytical re- 
sults) for distances between 1 nm and 10 /im for the 
two configurations are plotted in Fig. ^. Since typically 
T K, 10^^"^ s, the long distance regime {D ^ ct ^ 10 /im) 
will not be visible in these plots. For the plots, a Drude 
model is used with hojc = 5.9 meV and hup = 9.85 eV. As 
expected from the analytical results, the force in the in- 
plane configuration will be larger than that for the polar 
configuration for small enough distances {D < 10 nm). 
The discontinuity at D w 40 nm in the plot of the in- 
plane case depicts the change of sign. The analytical 
results are in pretty good agreement with the numerical 
calculations. 
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IV. THE HYBRID MODEL 

The Drude model is not very realistic. Although it 
describes rather well the diagonal part of the dielec- 
tric tensor (except of course for the effect of interband 
transitions, which are not very important here), the off- 
diagonal part of the dielectric tensor is poorly described. 
This is because the latter is dominated by interband tran- 
sitions. We therefore introduce a model (called "hybrid 
model" ) in which e^x is described by a plasma model. 
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FIG. 2: Absolute values of the magnetic Casimir force (per unit area) with the mirrors described by the hybrid model. Numerical 
results (solid curves) are compared with the analytical expressions (dashed curves) for the polar configuration (a), the in-plane 
configuration (b), the term AJ^f resulting from the longitudinal Kerr effect (c) and the term from the transversal Kerr effect 
at! (d). 



and where e^y is described by a single absorption line 
interband transition: 
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In real systems luq will be of the same order of magnitude 
as ujp. The off-diagonal element of the dielectric tensor 
at imaginary frequency can be obtained by the following 
Kramers- Kronig relation: 
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For this model, we will only have two different integration 
regimes; the long distance {D ^ c/ojp) and the short 
distance regime (D <^ c/cup). In the long distance regime, 
the integrals in Eqs. (jj) and (||) will be dominated by the 
range u < k±c <^ ujp. In this range, we can approximate 
the dielectric tensor by: 
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One then finds for the reflection coefficients: 
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With these approximations, we obtain for the magnetic 
Casimir energies and forces (per unit area): 
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As in the Drude model, the force in the polar config- 
uration will be negative. However, the total force for 
the in-plane configuration will be positive for the hybrid 
model in this distance regime. 

In the limit of short distances, one has to distinguish be- 
tween three different integration ranges while perform- 
ing the integrals in Eqs. M) and (P): region (i) where 
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by Eqs. (25^25^)!^ In regions (ii) and (iii), we will do 
the calculations without multiple reflections (i.e., we put 
rss(iw,ifc±) = rpp(iu;,ifc±) ^ 0). In region (ii), the di- 
electric tensor elements are given in Eqs. (pj), and the 
magneto-optical reflection coefficients are given by: 
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In region (iii), the dielectric tensor can be approximated 
by: 
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One then finds for the magneto-optical reflection coeffi- 
cients in region (iii): 
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With these equations by hand, we made an interpolation 
for the reflection coefficients in region (ii) and (iii). This 
makes the calculation of the integrals less labor-intensive 
because we can calculate the two regions at the same 
time. We have put: 
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With these expressions for the reflection coefficients, we gies and forces for the short distance regime. The result 
are finally ready to calculate the magnetic Casimir ener- 
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to define our integration range in more detail. In our 
calculations, the complete range of data extending from 
4 meV to 10 keV available in ||] was used, along with 
a Drude model below 4 meV as shown in Fig. The 
following parameters for the Drude model were found 
by extrapolation of the available data at low frequen- 
cies: fkOp = 3.5 eV and h/r = 19 meV. The quantity 
£2:2; (iw) — 1 calculated in this way is shown in Fig. |^ to 
decay roughly as uj~^/^ (for u) > h/r), so it can not be 
completely described by a Drude (or plasma) model. 

Experimental data for the ofF-diagonal element of 
the dielectric tensor is rather scarce. Some data for 
Re Sxyiuj) can be found in 0. They are shown in Fig. IJ. 
With the causality relation (|22|) it is then possible to cal- 
culate Exyiii-^)- Since we only have data available between 
0.1 eV and 6 eV, we had to perform the integration in 
Eq. ( p^ ) over this range. This, of course, is a rather rough 
approximation. The results of the calculation depicted in 
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for distances small enough. The force for the in-plane 
configuration is positive in this short distance regime too, 
so there will not be a change of sign for this model. 
We did the same numerical calculations as for the Drude 
model and compared the results with the analytical ex- 
pressions in the different regimes. The following param- 
eters were used in our hybrid model: Hujp = 9.85 eV, 
fkuo = 3.9 eV, — 1.5 10^^ and we have put lo* — 
2 exp(l)ti;p. In Fig. ^ the numerical and analytical re- 
sults for the absolute value of the magnetic Casimir force 
are shown to be in rather good agreement. 



V. NUMERICAL CALCULATIONS ON FE 

The Drude and hybrid model will not provide an ac- 
curate description for the dielectric tensor of the mirrors 
in a real system. This is because interband transitions 
will start playing a role at photon energies of a few eV, 
and these are not contained correctly in either one of the 
models. In order to obtain an estimate of the magni- 
tude of the magnetic Casimir force in such a real sys- 
tem, it is necessary to perform numerical calculations of 
Eqs. (^ and (H) where the reflection coefficients are cal- 
culated with experimental data for the dielectric tensor. 
In this section, we will present such calculations for a 
system with iron plates. Similar calculations for the non- 
magnetic Casimir force have already been performed for 
Al, Au and Cu § |. 

Experimental values for the imaginary part of exx{^) for 
Fe are given in ||^ . The diagonal element of the dielectric 
tensor at imaginary frequency can then be obtained by 
the causality relation: 
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I Lo' Im exx{<^') 



(32) 



Of course it is impossible to perform the numerical in- 
tegration over the entire interval [0, -|-cx)], so we have 



hybrid model (cf. Eqs. ( P4b| ) and (p8b|)). 

The magnetic Casimir force and energy are now calcu- 
lated by numerical integration of Eqs. (^) and (^. We are 
interested in plate separations between 1 nm and 10 /.tm. 
These separations correspond to frequencies in the range 
10~2 - 10^ eV, so we will have to perform the integra- 
tion between say 10~^ eV and 10* eV. Fig. ^ shows the 
resulting force and energy (per unit area) for the polar 
and in-plane configuration. In the short distance limit, 
the force decays as Z?"^ for the polar configuration and 
as I?~3 for the in-plane case. For long distances we find 
a power law for the polar configuration and D^^ for 
the situation with magnetization parallel to the plates. A 
change of sign of the interaction for the in-plane config- 
uration is also visible from the figure (the discontinuity 
at = 50 nm). The power laws differ (except for the 
in-plane configuration at short distances) from those ob- 
tained for the Drude and hybrid model. This can be ex- 
plained as due to the different behavior (because of inter- 
band transitions) of the dielectric tensor for Fe compared 
to that of the models. In view of future experimental in- 
vestigations of the effect, distances Z? > 10 nm are the 
most interesting. In this range, the effect will be great- 
est for the polar configuration. For two parallel plates 
of Fe (with infinite lateral extension) the force per unit 
area in this configuration is approximately 40 mN/m^ at 
D = 10 nm, and decays to 0.1 mN/m^ at = 100 nm. 
Whether such forces can be observed experimentally will 
be discussed in the next section. 



VI. EXPERIMENTAL SETUP 

Since it is hard to experimentally maintain two parallel 
plates uniformly separated by distances smaller than a 
micron, one of the plates is most often replaced by a 
lens-shaped mirror. Recently, a number of experiments 
has been performed using this geometry to measure the 
non-magnetic Casimir force with an AFM ||l^, |l^, p^ . 
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FIG. 3: The imaginary part of the diagonal element of the dielectric tensor evaluated at real frequencies (a), and the diagonal 
element as a function of imaginary frequency (b). 
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The real part of the off-diagonal element of the dielectric tensor evaluated at real frequencies (multiplied by lj^) (a) 
off-diagonal element as a function of imaginary frequency (b). 



For this plate-cylinder geometry, the Casimir force can 
be obtained from the plate-plate geometry by means of 
the force proximity theorem |T3] : 



AF = 27rRA£{D). 



(33) 



In this equation, R is the radius of curvature of the lens- 
shaped mirror, and A£{D) is the Casimir energy per unit 
area for the configuration with two plates. One has to be 
careful to distinguish between AF and AJ-; the former 
is the force for the plate-lens geometry, while the latter 
is a force per unit area for two parallel plates. With the 
numerical results from the previous section, we are able 
to estimate the magnitude of the magnetic Casimir force 
in the plate-lens geometry for Fe. If we take R = 100/xm 
and a distance D = 50 nm, a force |Ai^| « 10 fN is found 



for the polar configuration. In the in-plane configuration, 
the force will be two orders of magnitude smaller. Such 
small forces can probably not be measured with the AFM 
technique. However, sensitivities of 0.1 to 10 fN in "mag- 
netic resonant force microscopy (MRFM)" have been re- 
ported |l^. More recently, the detection of forces in 
the attoncwton range has been achieved ||l^, |l^ . 

A possible MRFM setup is already discussed in de- 
tail in Ref. A thin film 10 nm) of ferromagnetic 
material with hard magnetization is deposited on a sub- 
strate which is placed on a piezo-electric actuator. The 
lens-shaped mirror is attached to a cantilever by first 
depositing a small droplet of polymer on the cantilever, 
which can then be covered by evaporation with a thin 
10 nm) layer of soft ferromagnet (such as permal- 
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D (m) D (m) 

FIG. 5: Absolute values of the magnetic Casimir force per unit area (a) and the magnetic Casimir energy per unit area (b) 
between two iron plates (of infinite lateral extension). The solid curve corresponds to the in- plane configuration, while the 
dashed curve describes the polar configuration. 



loy). In this way, one is able to create the lens-shape, 
with a curvature radius of say 100 ^m. The distance 
between the samples can be controlled easily with the 
piezo-actuator. By applying an ac magnetic field, one is 
able to modulate the magnetization of the soft sample at 
the resonance frequency of the cantilever. This will gen- 
erate an oscillating magnetic Casimir force that causes 
the cantilever to vibrate. The deflection of the cantilever 
can then be measured with a laser. In this way, the mag- 
netic force (AF = Faf — Ffm) can be measured. The 
force resolution achievable using a freely vibrating can- 
tilever is fundamentally limited by intrinsic thermome- 
chanical noise. This force noise can be controlled by the 
geometry of the cantilever; one needs a high Q cantilever 
that is thin, narrow and long to obtain the best sensitiv- 
ity. With ultrathin silicon cantilevers, force resolutions 
in the attonewton range have been obtained jl6j. More 
information on the sensitivity of MRFM can be found 
in P, [16 . 



Since the non-magnetic Casimir effect is independent 
of the magnetization direction of the samples, only the 
magnetic contribution to the Casimir effect will be mea- 
sured by using this modulation technique. Parasitic elec- 
trostatic forces (caused by a difference in potential be- 
tween the magnetic samples) are also taken care of auto- 
matically in this way. The exchange interaction between 
the samples does not contribute much at the separations 
of interest (D > 10 nm). Another parasitic magneto- 
static interaction is the dipolc interaction between the 
ferromagnets. This dipole force can be made as small as 
needed by taking a ferromagnetic plate with sufflciently 
large lateral extension and sufficiently small thickness. 
The plate should also be as uniformly magnetized as pos- 
sible. With a plate of radius 1 cm, and a thickness of 
10 nm, this parasitic magnetostatic force can be esti- 



mated to be below 1 attoNewton. Interaction of the soft 
sample with the ac magnetic field will yield a signal at 
two times the modulation frequency, so this can be fil- 
tered out effectively by using a lock-in amplifier. Thus 
with MRFM, it should be possible to measure the mag- 
netic Casimir interaction without much influence of other 
effects. 



VII. CONCLUSION 

In this paper, the magnetic Casimir interaction discov- 
ered in Q was generalized to the case where the magne- 
tization is parallel to the plates. The calculations for the 
Drudc model in the short-distance limit were revised, and 
another model was introduced. The behavior of the in- 
teraction was discussed in the different distance regimes, 
and it is seen that the interaction in the two models de- 
cays quite different with interplate distance. Numerical 
calculations for a real system with iron plates were also 
presented. Here we used experimental data for the dielec- 
tric tensor of the mirrors. The results from this numerical 
work on Fe could not be fitted by one of the introduced 
models, because interband transitions play a prominent 
role in Fe, and these were not implemented correctly in 
the models. 

It was made acceptable that the new Casimir magnetic 
interaction can be measured with magnetic resonance 
force microscopy. However, to obtain an accurate com- 
parison of the theory with eventual experimental results, 
more work is needed on the theoretical side. A detailed 
analysis of the geometrical effects would be valuable; the 
proximity force theorem does not provide reliable esti- 
mations at a level of accuracy of a few percent. Also one 
has to consider other corrections already calculated for 
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the non-magnetic Casimir effect; e.g. surface roughness 
corrections would probably play an important role 
Finally, more experimental data on the off-diagonal el- 
ement of the dielectric tensor for several ferromagnetic 



materials is also necessary in order to obtain a better 
estimate of the magnitude of the interaction from the 
numerical procedure presented. 
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